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ABSTRACT 

Cosmological simulations consistently predict specific properties of dark mat- 
ter halos, but these have not yet led to a physical understanding that is generally 
accepted. This is especially true for the central regions of these structures. Re- 
cently two major themes have emerged. In one, the dark matter halo is primarily 
a result of the sequential accretion of primordial structure (ie 'Nature'); while 
in the other, dynamical relaxation (ie 'Nurture') dominates at least in the cen- 
tral regions. Some relaxation is however required in either mechanism. In this 
paper we accept the recently established scale-free sub-structure of halos as an 
essential part of both mechanisms. Consequently; a simple model for the central 
relaxation based on a self-similar cascade of tidal interactions, is contrasted with 
a model based on the accretion of adiabatically self-similar, primordial structure. 
We conclude that a weak form of this relaxation is present in the simulations, 
but that is normally described as the radial orbit instability. 

Subject headings: cosmology: theory — dark matter — halo formation 

1. INTRODUCTION 



In Henriksen (2006b, H06; 2007, H07) a theory of dark matter relaxation has been pro- 
posed that is based on a temporally convergent series solution for the Distribution Function 
(true phase-space density or DF). The relaxation is effected in a non-mechanistic fashion 
by maximizing the local Boltzmann function calculated from this DF. This procedure deter- 
mines all of the parameters in the DF and allows the density distribution p(r) and the pseudo 
phase-space density 0(r) = p{r)/a{r)^ to be calculated, as well as any other quantity such as 
specific angular momentum. Studying these results one finds that, starting from the position 
and corresponding slopes found in the simulations, either the density flattens rapidly or the 
pseudo density steepens rapidly within the next decade of smaller radius. There is at present 
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no convincing evidence for either of these trends in the simulations, although unfortunately 
this region is at or near the current resolution limit. This paper attempts to understand the 
proposed relaxation more intuitively by suggesting a simple mechanistic model. Such a view 
may be briefly summarized as 'nurture'. 

In Salvador-Sole et al. (2007, S07) a rather different explanation for the simulated 

structures is proposed. In effect one uses the Press- Schechter (1974) formalism as expressed 
by Lacy and Cole (1993) to calculate the instantaneous merger rate. Since most of the 
halo mass is added by many small objects, this is approximated by being smooth. This 
smoothness allows one to accrete the dark matter halo from post-recombination large scale 
structure in a sort of 'layer cake' fashion, where each layer may be deduced from the accretion 
flow at the epoch when it was added. The innermost layers correspond to the earliest times 
so that 'inside-out' growth is estabhshed. The predicted density profile compares well to 
NFW (Navarro, Prenk and White, 1997; Navarro et al. 2004) profiles for various masses over 
appropriate ranges and even better over the whole range of scales to an Einasto (Navarro 
et al., 2004) or Sersic profile (e.g. S07). The latter profile has a finite density at the centre 
but an infinite slope. Similar comparisons may be made for velocity dispersion and angular 
momentum (Gonzales-Casado et al., 2007). Such a view may be briefly summarized as 
'Nature', although there is also a kind of relaxation that leads to adiabatic self-similarity 
(e.g. Henriksen 2006b, H07) in the asumption of "smoothness". 

Unless it is also hidden in the asumption of 'smoothness', the latter picture leaves little 
room for 'thermodynamic relaxation' (ic maximum entropy) as an clement in the simula- 
tions. Nevertheless, whether or not such relaxation is present in reality in order to explain 
possible density cores remains an open question. Moreover the remarkable results of the 
'Via Lactea' simulation (Diemand, Kuhlen and Madau, 2006-DKM06; Madau, Diemand, 
and Kuhlen, 2008-MDK08) show elaborate sub-structure that is characterized by a definite 
mass spectrum. Such a hierarchy should be interacting tidally and by dynamical friction and 
so it provides a mechanism for relaxation (see e.g. the discussions in El-Zant et al., 2004, 
H07, H06 and Henriksen, 2006a). 

One indication of the possible presence of such relaxation has been provided by Hoff- 
man et al. (2007). They observe that the halo surface density at the NFW scale radius 
(that is psTs) remains constant during the evolution of an individual halo. This together 
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with virialization within allows them to conclude that (ps cc Vg . Some time ago, virial- 
ization together with strong dynamical inter-scale coupling were shown to be equivalent to 
constant surface density and virialization throughout a cascade of structures. This was in 
a completely different context (Henriksen and Turner, 1984; Henriksen, 1991 and references 
therein) where it was studied in the context of star formation. It was proposed that a hier- 
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archy of molecular clouds slowly evolved by coUisional and tidal interactions into the stellar 
Initial Mass Function (IMF) as the result of a kind of 'ballistic turbulence'. The observa- 
tion of Hoffman et al (2007) together with the power law substructure (DKM06) suggests a 
similar argument may apply to dark matter halos. 

In the stellar case it has been difficult to find any evidence for such an initial fragmented 
state, and so ironically the picture may be more relevant to dark halos than it is to luminous 
stars! However Hoffman et al. have not applied the idea of a dynamical cascade to all 
scales below r^. This will be the subject of the next section. In the third section we give 
a naive calculation that imitates the layer cake approach of (S07), by simply assuming 
adiabatic self-similarity as dictated by the primordial perturbation spectrum. Finally in the 
conclusions we discuss the likhhood of either sort of relaxation. We conclude that a weak 
form of cascade relaxation is present in the simulations that is not distinguishable from the 
radial orbit instabihty. We conclude further that the pure cascade evolution should lead to 
a sharp break in the density and pseudo- density trends as simulated currently in the next 
decade or so of resolution. Otherwise we are left only with the weak form of relaxation that 
is described in this paper as either 'nature' or adiabatic self-similarity, and which seems to 
be equivalent to the radial orbit instability. 



In this section wc apply the self-similar model of cascade structure introduced for molec- 
ular clouds in Hcnriksen and Turner (1984; HT84 subsequently) and summarized in Hen- 
riksen (1991; H91 subsequently). We describe first the 'initial' self-similar structure of the 
cascade and subsequently estimate its equilibrium radial structure using a local maximization 
of entropy (following Henriksen (2007)). 



2. 



Cascade Relaxation 



2.1. 



Initial Cascade 



We use a fractal description of the number of halo substructures according to 



n(£,r,) = (-) 



-D _ 



dN{e,rs) 
dlnie/rsY 



(1) 



Here n{i, Vg) is the number of substructures inside the NFW scale radius each having 
scale i, and is a cascade parameter that is usually described as the fractal index since 
it need not be integer. This is the true fractal view wherein the structures are separated 
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by discrete steps in scale which, in a self-similar fractal, are effected by a constant scale 
ratio. In fact equation ([1]) follows by recursion once the number of substructures of size £ 
in an object of size L is assumed to be n{l, L) = {£/L)~^. The transition to the continuum 
description is also indicated in equation ([T]). We set n(£, r^) equal to the derivative with 
respect to the logarithm of scale (so equal differences are equal ratios) of the cumulative 
number of substructures inside having scale less than i, namely N{i,rs). 

The description of the cascade used in this paper can appear confusing because of the 
intrinsic number of parameters and of several ad hoc notations that are used either for brevity 
or for correspondence with previous work. For the convenience of the reader we assemble 
here the various notations in table 1, even before their introduction in the text. There are 
really only two cascade parameters under our assumptions namely a and (3 in terms of which 
all other cascade quantities may be expressed. These are defined below just before equation 
(j3]) and by equation (fTTl) . We use generally a systematic ad hoc notation wherein stands 
for the negative power of the dependence on radius of the quantity x. The symbol D without 
a subscript corresponds to the fractal index of the cascade elements as above. 

Table 1: Definitions^ 



Parameter 


Definition 


D 


2 + (3 


Dp 


I -{a + (3) 




(5+a+/3) 
2 


Dp 


-(3 + a + /5) 


V 


3 + 2(« + /3) 


ii 


a 


K 


In {rs/Q 


L 


Cascade top scale 



°Ad hoc quantities in terms of Cascade parameters 

We define the scale of each substructure as a kind of gravitational correlation length 
(e.g HT84) over which the velocity dispersion a{t) and the total mass M{t) satisfy the virial 
condition 

- (2) 

We show below (equation fl27|) and subsequent remark) that a self-similar virial cascade 
of these sub-structures implies a hierarchical vector velocity that adds randomly to give a 
virialized structure on the scale r^, as it should. This local virialization assumption has been 
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found to hold recently in a general class of distribution functions (Wyn Evans and An, 2005, 
An and Wyn Evans, 2005). This is reasuring since the assumption has been used intuitively 
for some time in these cascade models, as noted in the introduction to this paper. 

In order to parameterize each substructure as composed of substructures plus a smooth 
distribution of particles on the scale L, we introduce the fraction f{i, L) as the fraction of 
mass M{L) contained in structures of the next sub-scale £. In a self-similar fractal it is 
natural to take this as f{i, L) = {i/L)°', where a is thus another parameter of the cascade. 
Consequently the mass of an object of scale L is given by 

By recursion this yields the mass of an object of scale i, in terms of the total mass inside rs, 
namely Ms, as 

M{e) = M,(-^)(^+"). (4) 
the mean density of an object of size i is therefore 

where ps = Ms/{4.-Kr^j3). 

In H91, the mean density is taken to vary as £~^p so that in the present notation 

Dp^3-{D + a), (6) 

which may be employed where convenient. Moreover a is thereby seen to be the 'aggregation 
index' ii as defined in H91. When this quantity equals zero, all of the mass at each scale 
is composed of identifiable substructures. We will refer to this condition as an 'aggregation 
cascade' for brevity. In H07 the index Dp = 2a. 

We will, for the purposes of describing the simulations with this model, frequently want 
to make the correspondence between scale £ and spherical radius r. In the spirit of the 
dynamical coupling between scales advocated in HT84 and in H91 we suppose that this 
correspondence follows by equating the dynamical time at r to the dynamical time of the 
maximum sub-structure of scale ir permitted at r, so that 

aH£rl ^ GM{r) 
£1 ' ^ ' 
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Hence recalling equation ([2]) for the scale f^, and using equation (jlj) for a sphere of radius r 
(rather than r^) to obtain M(r) = M(£r)(^r/'")~^^^"\ we conclude that 

(r/4)^''+"-'^ = 0(1), (8) 

and hence ~ This means that the maximum scale at any r is about equal to r itself. 

The first constraint on such a hierarchy from the simulations comes from the DMK06 and 
MDK08 papers where they present evidence for the substructure cumulative mass function 
N{Msub) in the form 

« (9) 



d In Msub 

where q is approximately in the range 0.9 < q < 1.0. In our model Msub = M{€) and by 
equation (jl]) d\n.M{t)/Ms = {D + a)d\n.i/rs so that by substituting into the continuum 
form of equation ([T]) we have our prediction for this result 



dN{M{l)) _ / 1 \ / M(£) \~^ 
dln(M(£)/M,) ~ V^ + aJ \ J 



(10) 



It is therefore already clear that such a model can only fit the simulations if a « 1. 
This would imply an approximate aggregation cascade (i.e. ii = a ~ 0) at least at the larger 
scales. This is not consistent with the 2 = simulation state (DMK06, MDK08), wherein 
they find that the sub-structure is only about five percent of the total mass inside the virial 
radius. However we must remember that at present we are describing the early state of the 
cascade where such a condition may have been more appropriate. More appropriate because 
although the mass function does not seem very sensitive to epoch in the simulations the 
mass fraction in clumped material increases with z (MDK08). 

The key to the subsequent evolution of such a cascade is the interscale dynamical cou- 
pling, which on our view promotes the relaxation of the system (HT84,H91) and is due im- 
plicitly to tidal interactions, which we take to include collisions with all impact parameters. 
This is quantified by comparing the dynamical time on a scale L, namely tdyn = L/a{L) to 
the collision time of the next internal subscale namely tcoii = l/{{n{i,L)/L^)TTi'^a{L)). A 
strongly coupled cascade would have tcoii/tdyn ~ 1, but we generalize slightly by introducing 
a third parameter f3 such that (vr is introduced for convenience) 

^=-(y)^ (11) 

SO that by using the various definitions 



or finally 
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Since this can only hold for a scale ratio unequal to 1 if the power on the left vanishes, 
our cascade parameters are again reduced to two (i.e. a and j3) as now we require 

D = 2 + P, (14) 

and hence 

Dp = l-{a + (3). (15) 

If (3 is negative the cascade becomes increasingly collision dominated at the small scales 
(see equation ( fTTi) ) and we expect the relaxation to be stronger there. Ultimately the sub- 
structure may be expected to merge into a continuum at small scales. 

The pseudo phase-space density (Taylor and Navarro, 2001) may also be calculated in 
the cascade as p{i)/a{i)^ from equations ([5], El Hj) to be 

= 0.(-)-(^+"+^)/^. (16) 

Here 0s = Ps/c^l and recall that we may expect to replace ^ with r on average. For brevity 
we may refer to (5 + a + /?)/2 as simply D^. 

In the strongly coupled, aggregation cascade (/3 = 0, a = 0), this varies like r~^-^ while 
p oc (ie constant surface density, HT84). These relations should hold everywhere in the 
un-evolved cascade, not just at any particular radius such as the NFW scale radius (Hoffman 
et al 2007). 

However this simple description does not fit the 2 = state of the simulations. This state 
requires that at a radius somewhat smaller than the NFW scale radius we have approximately 
a + (3 ~ —1, so that Dp ^ 2 and 2. In addition a ^ in order to fit the sub-structure 

mass function . As the simulations are continued to smaller radii (down by about two orders 
of magnitude) Dp decreases to about 1 and slowly thereafter (Navarro et al., 2004), 
appears to remain close to 2 (Austin et al. 2005; Dehnen and McLaughlin, 2005) and a 
remains close to zero to the resolution limit (DKM06). 

We do not expect that relaxation by cascade is relevant much beyond the NFW scale 
radius where the density profile is steeper that and the 'nature' of the environment 
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probably predominates (see e.g. below and Henriksen and Widrow, 1999). In contrast, 
the evolution toward smaller scales should indeed manifest relaxation by this clump-clump 
mechanism. We thus turn next to examining whether or not this simple cascade model can 
evolve to smaller radii as expected from the simulations. 

Such an examination was carried out also in H07, but the calculations were rather 
formal. They were based on a renormalized, 'adiabatically self-similar', distribution function 
together with the maximization of the local entropy density. It is this assumed maximization 
that takes the place of the detailed dynamics of the tidal interactions. In the next section we 
give a simpler version of this discussion based on the preceding cascade structure. Similar 
conclusions to those of H07 are reached nevertheless. 



2.2. Cascade Evolution 

The actual dynamical implementation of a cascade of tidal interactions or collisions 
escapes a precise description for the moment. To avoid this impasse, we use the idea that 
these processes move the cascade through a succession of local equilibria (e.g. H07). Each 
equilibrium on a scale £ is taken to be associated with a maximum of the Boltzman H 
function. In Henriksen (2007) this was discussed in terms of parameterized distribution 
functions, whose parameters varied with scale. We attempt here to give a simpler version of 
the model in terms of the sub-structure cascade detected in the simulations. 

As a result of this assumption the details of the interactions are ignored. They enter im- 
plicitly however through equation (ITT]) . Should /?>>!, then tidal collisions and dynamical 
friction would not play a role in the dynamical evolution below r^. However our evolution 
below gives the range of a + (3 and together with a small this range requires to never be 
far from unity in order to approximate the simulations. And in fact we shall see that the 
relaxation detected in the simulations requires j3 < 0. 

In Henriksen (H91) it was suggested that the cascade evolved on a given scale when a 
dimensionless time r(t, £) reached a critical value tj on that scale. In the stellar case this 
heralded the transformation from molecular cloud sub-structure to a pre-stellar 'core'. How- 
ever this depended essentially on dissipation, which is not present in the present application. 
In this context we expect the parameters of the initial cascade to be modified on each scale at 
some dimensionless time due to tidal transfer of energy and angular momentum. We expect 
this process to go faster on the smaller scales and to transfer energy, angular momentum 
and mass to larger scales. 

As a simple implementation of this last remark, in order to provide some expectation 
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of a timescale for the equilibrium calculation to follow, we may estimate the form of tj by 
dimensional analysis of the initial cascade. Given a cascade characterized only by power 
laws in radius for both the density and the pseudo-density, their dimensional factors (say 
A, n respectively) may be taken as the dimensional constants that define the cascade self- 
similarity. These may be combined with the gravitational constant to yield two expressions 
for a local time scale, namely from the density 

from the pseudo-density as = = l/(G0sr^)(£/r>,)~*^^~^<*)i(:. The dimensional 

constants have been written here implicitly in terms of the values of density and pseudo- 
density at Ts, namely A = p^rf'' = p{i)i^'' and /i = (psf^''' = 4'{i)i^'^. In order to obtain 
a single estimate of the evolutionary time scale we would like these two estimates to agree. 
Happily this requires = 3 — Dp/2, which is in fact identically satisfied by their definitions 
in terms of cascade parameters a and f3. They also agree numerically under our conditions. 

Other possible time scales for the cascade, such as the global dimensionless dynamical 
time {o's/f^s)t, are not independent of these dimensional factors. This latter dynamical time 
is in fact just y/Uplt, which is Tp above at the scale r^. 

Setting the common r defined above equal to a value tj independent of scale estimates 
the initial evolution of the cascade. We see that the scale of the evolved cascade increases 
with a positive power of t, the age of the system. Hence, descending in scale, we expect to 
find an increasingly evolved cascade at a fixed t. The numerical value of tj may be estimated 
to be of order unity since it is numerically equal to y/Gpst and we expect the relaxation to 
occur in several dynamical times on the scale r^. Thus the time (from the origin of the 
cascade) to significant evolution on a scale i, say tj, is given by 

tAi) = -^{-r^''. (17) 

Hence as either we descend the sub-structure scale, or as Dp decreases in time (see below), 
the relaxed state becomes less dependent on scale. Such a picture agrees with equation ffTTl) 
when the initial (3 is negative. 

The preceding estimate can only be approximate as it does not take into account the 
adiabatic evolution of the cascade parameters. To achieve this, as remarked above, the 
actual dynamics of this evolution is replaced by a statistical argument based on maximizing 
the local entropy density (Boltzmann H function: see e.g. H07 and discussion above). By 
applying it locally in radius we hope to find the radial dependence of the evolution in the 
cascade parameters a, (3. In this connection we should recall the work of Hansen et al. (2006) 
who studied the velocity distribution in dark matter simulations, as well as similar work by 
Merrall and Henriksen (2003). In these papers evidence for relaxation (near Gaussianity 
and universality) is found in colliding and collapsing systems. In the Merrall and Henriksen 
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paper, the radial velocity distribution function is fitted to a Gaussian (with distortions) 
under many conditions, while in Hansen et al both the radial and tangential velocities are 
fitted to a Tsallis entropy distribution. Some anisotropy is detected in the latter paper 
that argues against complete thermodynamic relaxation, but relaxation is still present. This 
anisotropy is liable to be dependent on net angular momentum in the system which tends to 
slow the approach to equilibrium. We ignore this effect in the nurture calculation, but not 
in the nature calculation, where the relaxation is indeed slower. 

To proceed we estimate the phase-space volume fl available to substructures of scale £ 
inside the sphere of radius r as (for distinguishable particles) 

mr)^{ 3 )"^^'-^ (18) 

Here n{£,r) is the number of sub-structures on scale £ inside a sphere of radius r, acm{£) is 
the average velocity dispersion of the corresponding element relative to the centre of mass of 
the halo and is a convenient element of phase space with which to scale the phase space 
volume. 

It is in fact appropriate to take as a volume close to the largest scale in the cascade, 
Vs, since the smallest scale is not well determined. This leads to the inverted form for the 
phase-space volume above in order to keep the entropy positive. Thus we write h = £o(To 
where £0 and cTq are quantities on a scale close to the top of the cascade. Using equations 
(I2I Hj) and the various definitions 

-W = c..(-)^, (19) 

rs 

so that we obtain h by applying this last result to £0 and then multiplying by £0 as 

/ £n N 3 + a + /3 , . 

h = rsas{-)^. 20 

rs 

We expect £0 to be somewhat smaller than to form the top of the cascade. The advantage 
to this choice for h over simply the product of and ag is that it allows for the evolution 
of the cascade with varying parameters a and jS. We find below that h increases during the 
evolution so that the cascade indeed grow from small scales to large scales. 

We calculate acm{£) as the result of A^(^) steps (not to be confused with the cumulative 
mass function) of a vector random walk between scales £ and r^, each step of size TZ. We 
should first note that equation (jlj) can be written in terms of the N*^ step as 



M{N) = M,7^^(^+"^ 



(21) 
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where N is related to i by 

N(0 - (22) 
Then we write the halo centre of mass velocity on the scale namely as 

vn = vn~i + Awat, (23) 

where Awtv is the randomly oriented velocity step from — 1 to A^. We take this once again 
to have amplitude equal to the virial velocity on that scale, namely \J GM{N) /TZts- Hence 
by squaring and averaging over all directions so that we may drop the cross term we obtain 

<yl = <^l-i + (24) 
and so by recursion to as the largest scale where = 

<yUN) = al (l + Sjlf 7^^■(^+"-l)) . (25) 

We have set GM^/r^ = cr^. 

On summing the preceding geometric series we find finally 

a,M^ = a'lF\n,a,(3), (26) 
where we use the ad hoc definition 

F2(7^,a,/3) ^ (^l+7^(^+°+^)( ^ )J . (27) 

Note that on setting N{r-s) = we obtain F = 1 and thus equation (J26i) merely restates the 
virialization on the scale r^. 

We will normally choose 7^ = 1/3 arbitrarily in our examples, since changing this value 
merely changes the scale associated with a given A^. 

The mean population n(£, r) is estimated as n{i)P{£, r) where P{i, r) is the probability 
that a structure of scale i should be found inside a sphere of radius r. We write this as 

P(£,r)= ^ ^^CM gFZE!^ ^28) 

Zj 

where the 'partition function' Z is given approximately by the integral 

Z(7^,a,/5)= dre 2gm(.)- 

= \^^ duu-'^e-^'^. (29) 
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We use the ad hoc notation 



c-'^^^n^. (30) 



2 

and u is a. convenient combined cascade parameter defined as 

u = 3 + 2{a + p). (31) 

The convenient variable u = {t/ts)"^ . 

There is a degeneracy between the parameters a and (3 as they only appear below as 
a sum or equivalently as {v — 3)/2. The partition function can be evaluated in terms of 
Whittaker functions, but we shall not need this in the first approximation. 



r as 



We now construct a local 'entropy' for the gas of sub-structures d in the sphere of radius 

^(£,r) = A;lnfi(£,r), (32) 



and by assuming that this attains a maximum as a function of the cascade parameter v and 
X = r /ts for a fixed scale (equivalently A^) we obtain the equation (cf H07) 

dx dJnS' ^ ' 

which computes the evolution t/(a;; N). Explicitly this takes the form 

du _ Se-""" + C^uA* 

dlnw ~ (3/2)e-""'(9^1nF2 + K/2) - {C^w + d^C^ + e-'^'^iflnTZ + In Z)) A*' ^ ' 

In this equation we use the previous 'ad hoc' notation plus another as 

A*{x;R,u) = -3w + ^lnF^ + ^K{u + 3), (35) 



where w = —\iax, and = In {rs/C-o)- 

Equation (1341) is a rather complicated differential equation for z/(x; N) in which a singular 
point (0/0) may appear. Fortunately, except possibly at the singular point, we can simplify 
it considerably by realizing that the derivative of the logarithm of Z is small and that is 
never very far from 1 for u > 1. Adopting = 1 also simplifies substantially. This gives 
us a working equation in the form flMl) but without the terms in InZ and InF^. Moreover 
from equations flSUj 1271) we have 

C2 Pi ^exp(^^^A^ln7^). (36) 
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An examination of the singular point in the simphfied equations shows that in general 
it is a focal point and solutions can not be expected to pass smoothly through it. It seems 
unlikely that adding the small terms omitted above will change this behaviour. 

It transpires that the initial condition that we wish to apply to equation (IMI) is z/ = 1 
at w = 0. This ensures that near we have a + (3 = —1 so that Dp = = 2 as simulated. 
Since we expect a to be small, we can assume that /3 ~ — 1. This means that the initial 
cascade is 'uncoupled' at larger scales (cf equation [TTl) . 

Equation flMl) allows us in principle to find how each scale [N or i) evolves to smaller 
r, as the local entropy attains its maximum. Should u increase towards 5ora + /5 = +l, 
then a density core has formed {Dp = 0) at that scale and the corresponding radius. If 
this evolution occurs before encountering the singular point, then our approximations are 
justified. Moreover because of the focal nature of the singular point, there is unlikely to be 
a unique critical curve that gives the 'physical' behaviour. 

We see in figure ([1]) that u does indeed evolve in this fashion over a range of scales. 
As was found in H07, this behaviour also implies that takes on the value 3 as the core 
flattens. 

More particularly in figure ([1]) we see that u for all scales remains close to its initial value 
while w increases by unity, after which it rapidly converges to universal linear behaviour in 
w = — Inr/rg. This behaviour requires the negative of the density slope Dp to pass through 
the value 1.5 in the same interval, after which it declines linearly with w for all scales. The 
negative pseudo-density slope D^ rises slowly over the same interval (more rapidly for the 
top of the cascade) after which it increases linearly with w for all scales. This logarithmic 
behaviour of Dp and D^ was also found in HOT and was shown there to be a reasonable fit 
to the empirical results (Navarro et al, 2004)over a limited range in r. Ultimately at smaller 
r it appears to vary too rapidly relative to the simulation. 

We note that this behaviour z/(a;) is all on one side of the singular point of the differential 
equation, but such solutions yield the physical regime for our parameters. In fact only for the 
case N = 2 (see e.g. figure does the focal point fall into the interesting region 1 < z/ < 5, 
when it occurs at w ~ 1.8335, u ~ 2.2273 (where the slope of the curve Dp becomes infinite). 
The variation in z/, and hence in Dp and D^, can be slowed by reducing the upper scale of 
the cascade (increasing K), but then the singular 'region' may be reached before z/ = 5 
as for N = 2 above. This behaviour is indicated in figure ([2]). We see that as the top 
of the cascade is reduced in size , the smaller scales (larger A^) slow down their evolution 
substantially. Indeed it would be difficult to distinguish them from the nature case for 
n = 4.5 shown in the lower right panel of figure([T]) over a limited range in scale. 
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Fig. 1. — The left panel shows the evolution of v with log of the reciprocal radius for three 
different scales. The top line is for iV = 4, while the bottom line gives the behaviour both 
for = 8 and N = 2. The right panel shows the corresponding density slope Dp with 
N = 8, 4, 2 respectively from the top curve to the bottom curve. The bottom left panel 
shows with = 2, 4, 8 respectively from the top curve to the bottom curve. All cases 
have 7?. = 1/3 and K = 1. The bottom right panel gives Dp from the 'nature' calculation 
below for 6 = 5 (dashed line) and 6 = 4.5 solid line. 
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Fig. 2. — This figure shows the same cases as found in the upper right panel of figure ([T]) 
but with the value K = 2. One sees that for N = 2 (bottom curve), which is now at the top 
of the cascade, the focal point behaviour is encountered. The curves for = 4 and = 8 
(middle and top curves) are substantially slowed relative to the K = 1 behaviour. 



Although the preceding suggests that the interacting cascade starts well below in 
order to agree with the simulations, the qualitative agreement between this simple model 
for cascade relaxation and the predictions of the renormalized series for the distribution 
function ( H07) is marked. We have simply added here an explicit mechanism by which the 
local entropy may maximize, inspired by the structure cascade found in DKM06. 

Our conclusions in this section are moreover the same as in H07. If relaxation can 
occur in such a fashion as to maximize the local entropy, then there should be a flattening 
density core and a steepening pseudo-density in relaxed dark matter halos. Moreover this 
should happen relatively rapidly for scales near the top of the cascade, although this upper 
limit to the relaxed cascade may be on a substantially smaller scale than Vg since it develops 
first at small scales. This allows the model to agree with the current simulations over a 
limited range. We observe further that the likely evolutionary path indicated by figure ([2]) 
is to descend the largest scale curve (A^ = 2) until the infinite slope in the density index is 
reached {w ^ 2). At this stage the radial scale has been reduced so that the next curve to the 
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right in the figure becomes relevant by a horizontal displacement. This continues until the 
universal behaviour is encountered at small scales, which is rapid compared to the 'nature' 
results also shown in figure ([T]). 

In the next section we turn to an analytic approximation to the 'nature' argument based 
on adiabatic self-similarity (H06) for the relaxation and smooth accretion. 



3. Dark Matter Halos from Cosmological Accretion 

In Manrique et al. (2003), Zhao et al. (2003) and Salvador-Sole et al. (2007; S07) the 
importance for the structure of dark matter halos of smooth accretion from cosmological 
conditions was emphasized. Indeed it was shown, particularly in the latter paper, that the 
assumption of smooth 'inside-out' accretion based on the accretion of small objects from 
the Press-Schecter paradigm (Lacey and Cole, 1993) allowed the simulated structure to be 
reproduced (i.e. NFW or Sersic density profile, Ms{rs) correlation) over an impressive range 
of scales. 

The key idea in these studies is that the density profile adjusts to current accretion con- 
ditions in a sequential 'layer-cake' fashion where the inner regions are fossil records of early 
cosmological times ('inside out' accretion). Moreover, whenever there is a reasonable approx- 
imation to a power-law spectrum for the cosmological perturbations of the form P{k) oc A;", 
S07 found that the usual models for self-similar relaxation (Fillmore and Goldreich, 1984; 
Hoffman and Shaham, 1985; Henriksen and Widrow, 1999; Le Delliou and Henriksen , 2003) 
do convert the cosmological accretion (Lacey and Cole, ibid) into the simulated density 
according to (our notation, with n an index independent of the previous section) 



o..2.-^-^. ,3r) 



Here 6 = 4 around a local maximum (Hoffman and Shaham, 1985) and 6 = 5 around an na 
peak (Henriksen H06). 

Thus even in this picture that relies directly on cosmological 'nature', there is some 
relaxation occurring. This has been described as "relaxation in an accretion bath" by 
Salvador-Sole (private communication) but in the language of Henriksen (H06) this may 
be considered as a kind of 'adiabatic self-similarity'. That is, the global constants deter- 
mining the self-similarity are dynamically changing. In the original picture (H06) this is 
due to local dark-matter relaxation, which we described in H07 and in the previous section 
by maximizing the local entropy. Here it is due rather to the changing cosmological infall 
with scale (and hence cosmological time) that dictates the local self-similarity through 
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equation fl37j) . Since the self-similarity is itself a result of relaxation, this is really a mixture 
of 'nature' and 'nurture' as remarked above. 

Just as in the previous discussion of Cascade relaxation, we include a table that defines 
the central quantities of this section. 

Table 2: Definitions^^ 



Quantity 



Definition 



P{k) 
n{k) 
k 



Q/Q* 

B,A 
a{k) 




"Adiabatic Nature quantities 
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Now in fact the usual approximation to a cold dark matter (CDM) power spectrum can 
not be universally described by a power-law, even of varying index. We have for example 
from Peebles (1993, p626) that 

{Q/A)Ak , , 

where 



A ^ m^^PC ^ IbAMpc 

4.7 

and Q is a scaling constant. If for convenience we choose to represent this as 



j^Mpc^ ^l7.^Mpc\ (39) 



P{k) = ^{Akr, (40) 

so that 

2 In (i + Afc + EF)- In (g/g-) 

n{k) = 1 , (41) 

then we see that one can only have a smooth adiabatic variation in n if one sets 

In g/g* = 2 In (2 + 5/^2). (42) 

This ensures a smooth passage through kA = 1. We observe that as ^ oo we have n —>■ —3 
slowly, while it goes slowly to = 1 as — * 0. Thus equation (1371) already ensures a flattened 
core and a power law density Dp = 12/(5 + 1) at infinity. 

But equation ( !37l) also allows the density profile to be calculated in detail as a function 
of A; ~ 27r/r, and consequently the variation with radius. Moreover (f) ^ p / {^J G pr'^)'^ gives 
D^{r) = — Dp(r)/2 + 3 as in the preceding section, so that this is also calculable. Moreover 
appealing to the general self-similarity relations found for example in HOT, Henriksen (H06, 
2006b) and in Henriken and Widrow (1999; S there should be read as 1/a = 2/ Dp) we have 
for specific angular momentum, mass and radius respectively 

M = M{X){^f-^'', (43) 
X 



(at)2/^P ' 



Here X, Z (not related to the partition function of the previous section) and A4{X) should 
be regarded as constants. Strictly speaking they are Lagrangian labels for different regions 
of the halo (just as are comparable factors in the density and pseudo-density), but in the 
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spirit of adiabatic self-similarity used as a description of a self-similar cascade (see e.g HT84) 
they may also be regarded as constant between the scales. All of these variations may be 
found based simply on equation fl?r|) . 

We therefore define an adiabatically varying index for as 

Dp{k) = -A + Dp{k). (44) 

We see moreover that near where Dp ^ 2, it is clear that Mg cx and that cx t, that 
is to the age of the system. Similar relations are found in simulations (e.g. S07; Zhao et al. 
2003; MacMillan, 2006; MacMillan, Widrow and Henriksen, 2006). 

This approach yields an analytic construction of a dark matter halo, based ultimately 
on equation fl57|) . It develops that the calculation with 6 = 4.5 gives a better fit to the NFW 
halos than either the case 5 = 5 or 5 = 4 so that we only show results for that case in figure 
([3]). When 6 = 4 the density increases too rapidly compared to the various NFW variations, 
while for 5 = 5 it is too slow. 

On the top right panel of figure ([3]) we see how well this approach can be made to fit the 
empirical profile of Navarro et al. (2004) (the middle curve) while the upper dashed curve 
is the original NFW profile, which is too steep. The lower curve has been calculated using 
equations fl37|) and fj4T|) . This curve agrees in essence with S07 over this limited range. We 
see moreover that increases slowly towards the value 3 once again while Dj2 increases 
slowly towards 4 as increases, both varying as Dp with r. The mean rms tangential 
velocity thus becomes linear with radius in the core and rolls over to ~ r°'^^ near r^. There 
is a remarkable similarity in the predicted evolution for D^ and Dj2. 

The interesting result is found in the comparison with figure (JT]). In that figure w = 
— In X is basically indentical with In k above. We see that there is qualitative agreement 
between the two approaches, but that the variations below occur much more rapidly 
in the maximum entropy approach. They occur in fact too rapidly to fit the numerical 
simulations over any substantial range (cf H07) if the relaxation begins close to rg. This is 
less pronounced if the relaxation happens mainly at smaller scales. The divergence is best 
seen in the bottom right panel of figure ([T]). This shows the variation in Dp for the present 
'nature' calculation, while the upper right panel shows it for the 'nurture' calculation for a 
cascade beginning at e'^Vg. Figure ([2]) shows the results when the cascade begins at e'^r^. 
All variations are ultimately linear in the logarithm of scale, but the slopes are very different. 

We must therefore conclude that this 'nature' mixed with the adiabatic self-similar 
relaxation describes the simulated halos very well, at least over a limited range in scales. 
Maximizing the local entropy in the 'nurture' approach can also be acceptable over a limited 
range in scale if the relaxation begins well below r^. However it predicts a rapid ultimate 
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Fig. 3. — The top left panel shows n{k) as found from equation ( HTi) . The top right panel 
shows a comparison between p{k) found from n{k) and equation (1571) and the NFW profiles, 
all normalized at k = 1. The middle curve is from Navarro et al. (2004), the upper curve 
is the original NFW profile and the bottom curve plots p from the present calculation. 
The bottom left panel gives the index D^{k) while the bottom right shows \Dj2{k)\. The 
calculations are presented for 5 = 4.5. 

divergence from the simulated behaviour. There is no evidence for such 'thermal relaxation' 
in the simulations at present. 
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4. Discussion and Conclusions 

The preceding sections allow us to conclude that the simulated dark matter halos are, 
up to the presently available resolution, consistent with the primarily 'nature' (S07) adia- 
batically self-similar smooth growth from the CDM spectrum (e.g. [3]). In our semi-analytic 
treatment of the previous section, we imitate this growth using the notion of adiabatic self- 
similarity (Henriksen, 2006b) to connect it to the CDM spectrum. In S07 however this 
was done using the assumption of smooth inside-out or 'layer-cake' growth according to the 
Lacey and Cole (1993) prescription. This employs a Press- Schecter (1974) type argument for 
the CDM spectrum, together with the statistical treatment of Bond et al. (1991). We con- 
clude that this may be done more simply using the notion of adiabatic self-similarity, which 
does involve a measure of relaxation. Moreover, although it happens slowly, this approach 
predicts a central core [n — > —3) rather than a cusp. 

What we have termed variously 'nurture', 'thermodynamic' or 'maximum entropy' relax- 
ation based on a cascade of interacting structure (H07 and above) agrees on the qualitative 
trends with the 'nature' view . However it evolves at small scales towards a fiat density 
and a steep pseudo- density more rapidly than is found either in the 'nature' discussion or 
in the simulations if it is assumed to extend close to ([1]). This discrepancy is less pro- 
nounced when the relaxation extends only to smaller scales, however (e.g. figure([2])). If we 
believe that cascade relaxation is playing a role in the halo evolution, then equation fill I) 
and the initial /5 ~ — 1 tell us that the relaxation should indeed be most effective at small 
scales initially. In time, as f3 becomes positive, the relaxation should move to large scales in 
agreement with equation ffTTl) . 

The only relaxation that is visible in the simulations at present is in the density profile 
within two decades or so of r^. The phase space pseudo-density shows no relaxation over 
this same range. This behaviour can be fitted over this limited range by either of the above 
relaxation mechanisms, but both mechanisms predict stronger small scale relaxation for 
which there is as yet no evidence in the simulations. 

One wonders whether the adiabatic self-similar relaxation that seems to fit the current 
simulations well, is in fact due to a form of cascade relaxation. Recently a strong case has 
been made (MacMillan, Widrow and Henriksen, 2006 (MWH06) and references therein) that 
the relaxation in this region near is due to the Radial Orbit Instability (ROI), so we should 
reconsider this in the present context. 

The onset of this instability was found in the MWH06 paper to coincide with the 
development of the mean square specific angular momentum at r into the Keplerian form 
oc GM{r)r. Let us suppose that the angular momentum perpendicular to the radial direction 
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takes the form £^0"^ where we use the same notation as in equation ([7]). Then the onset 
of the ROI is marked by the condition £^cr^(£r) ~ GM{r)r which on taking ir ^ r becomes 
essentially the cascade coupling condition of equation (^^. The explanation offered in the 
MWH06 paper is in terms of a resonant interaction between a 'bar-like' density perturbation 
and 'particle' (perhaps a sub-structure) orbits. In the present context we may see this as the 
largest asymmetric sub-structure at r. We expect it to be composed of sub-structures and 
to be interacting with such objects in the environment. 

It s interesting to note incidentally that this interaction is a kind of coupling between 
the radial infall and the transverse orbital motion, which delivers free-fall energy to the top 
of the cascade. The Cascade is really gravitationally driven turbulence. Ome might say 
that the mode k ^ 1/r parallel to the infall is coupling to the transverse mode kr ~ 
perpendicular to the radial infall to produce this turbulence. 

Of course these are all speculative, order of magnitude, arguments. Nevertheless we are 
inclined to suggest that the simulated halos, hence also the adiabatic self- similarity relaxation 
that works well for the current simulations, is due to the ROI. This is in turn a mechanism 
for the conversion of radial infall into transverse orbital motion at r. This transverse mode 
represents the top of the cascade at each radius. We thus identify the ROI with a weak form 
of cascade relaxation at large scales. This weakness is in accord with equations ffTTl) and f|T7j) 
that both predict (/3 < in the first case) stronger relaxation at small scales. The stronger 
cascade is expected to be at scales below the resolution of the current simulations. 

We therefore conclude that cascade evolution may be at work both in reality and in the 
simulations. We conclude tentatively that 'Nature' type adiabatic self-similarity is the weak 
form of the cascade relaxation and that it is equivalent to the ROI. Our predictions that 
would confirm the presence of cascade evolution are that there is a halo density core rather 
than a cusp and that the pseudo-density power law should break at higher resolution in both 
scale and mass. In the Via Lactea run this region is just inside their reported convergence 
radius at about Ikpc. There may be weak evidence in figure 1 of DKM06 that a density 
flattening is not excluded. Should this density and pseudo-density behaviour not appear 
down to say lOOpc in a Milky Way type halo, then the strong cascade relaxation is not 
present in the simulations. It is then probably not present in reality, unless for some reason 
the necessary interactions at a distance between sub-structures is discriminated against in 
the simulations. DKM06 do observe that at lower resolution than that of the Via Lactea 
run, the small sub-halos are poorly resolved. 

This work was supported in part by the Canadian Natural Science and Engineering 
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